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Abstract 

, We show that the one loop chiral corrections for heavy-light mesons in potential model can 

I explain the small mass of 1)5(2317) as well as the small mass gap between Z)s(2317) and £'(2308) 
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The recently observed 1^3(2317) 

HQS 

, which is a very narrow resonance (F < 10 
MeV) decaying into D^tt^, is thought to be the missing bound state with = 0^ of the 
heavy-hght system. This picture of i5s(2317) composed of a heavy quark and a hght valence 
quark fits well with the heavy-quark, chiral symmetries that predict parity doubling states 
(0~, 1~) and (0+, 1+), with the interparity mass splittings in the chiral limit given by the 
Goldberger-Treiman relation y, ^1. The subsequent observation of l"*" state D»(2460) 

nnn 

yj, 12, |3|| strongly supports this picture. 

On the other hand, the two-quark picture of the resonances does not play well with the 
potential model calculations, which generally predict substantially larger mass and width. 
According to the potential model calculation in Ref. I?! the mass and width of Ds{0^) are, 
respectively, 2487 MeV and a few 100 MeVs, with the width depending on the light-quark 
axial coupling. While the narrow decay width can be understood by the observed mass being 
below the threshold of the strong decay channel DK and the isospin symmetry breaking, 
the substantially small observed mass is puzzling. 

Furthermore, this anomaly in the observed mass became more peculiar when the Belle 
collaboration observed [8J non-strange 0"*" state D(2308), whose mass is surprisingly close to 
Z^s(2317). The potential model predicts the mass splitting between these states to be 110 
MeV. These peculiarities in the observed masses led to many models for the new resonances. 



including, for example, four-quark model |9|,ll0[, DK molecule models 



ll|, and unitarized 

meson model [1^. It is thus very important to clarify the nature of the newly discovered 
resonances. 

The quoted numbers of the potential model calculation are based on Coulombic vector 
potential and a linear scalar potential. Modifications of the employed potentials might 
remove the anomaly, but Cahn and Jackson ^| showed that, as far as the vector potential 
is kept Coulombic, it is unlikely that the observed decay width and mass pattern of the 
resonances can be obtained from a potential model. 

This suggests that the potential model be missing an essential physics of the heavy-light 
system. Indeed, the conventional potential model does not sufficiently take into account the 
chiral symmetry breaking nature of the QCD vacuum, with the chiral symmetry breaking 
encoded only in the light-quark constituent masses of the model. Since the light valence 
quark is chirally active the heavy-light mesons can couple to the quantum fiuctuations of 
the Goldstone bosons of the QCD vacuum. This suggests that potential models must be 



2 



augmented by chiral radiative corrections. 

In this paper we calculate chiral radiative corrections for the bound state energies of the 
potential model, paying particular attention to the mass splittings of the parity doubling 
states. Our main result is that chiral corrections are large, comparable at least to l/M^ 
corrections in D mesons, where denotes the charm mass, rendering their inclusion to 
the potential model mandatory. Furthermore, for the parity doubling states, they tend 
to narrow the interparity mass gaps, and this effect is stronger in strange system than in 
nonstrange system, with a robust prediction of the mass Gap = [m(D(0~*")) — m(Z}(0^))] — 
[m(L>,(0+)) - m(L),(0"))] ^ 90 MeV (at the axial coupling qa = 0.82) that is consistent 
with experiment. 

The potential model of heavy-light system j^] is based on the chiral quark model Q| , 
with the Lagrangian reading 

£ = ^\ido - H)<if (1) 

with \1' = (m, d, s) denoting the light quark fields and the Hamiltonian given by 

H = Ho + ^H, + ... (2) 

where M denotes the heavy quark mass. The leading Hamiltonian Hq in the heavy quark 
mass expansion reads 

^0 = 7°(-^ y + m) + V{r) (3) 

with the potential given in the form 

V^(r) = M + 7°V;(r) + K(r), (4) 

where Vg and Vy denote the scalar and vector potentials, respectively, and m = rriiSij denotes 
the constituent quark masses. The energy spectra of resonances are obtained by solving the 
Dirac equation of Hq, followed by time-independent perturbations of the sub leading terms. 
The free parameters of the model are fixed by a global fitting of the predicted masses to 
those of the observed resonances. 

In this framework the chiral symmetry breaking of QCD is encoded only in the con- 
stituent masses of the light quarks, and we shall see that this is not sufficient enough. This 
inadequacy of the model can be easily remedied by noting that the effective Hamiltonian is 



based on a truncated chiral quark model. In chiral quark model the light-quark-Goldstone 
boson interactions are described by an infinite tower of derivative expansions, but the term 
responsible for the one-loop corrections is the following axial coupling 

= ^^^.Yl,^jd,U,j + 0{U') (5) 
where qa is an axial coupling constant, and 

A, = l{ed,^-^d,e) (6) 

with ^ = e*"/2/.^ where U = Yll=i tt^A", A'^ the Gell-Mann matrices, and = 93 MeV. 

We note that the inclusion of the axial term (jS)) in the potential model Hamiltonian should 
not be unexpected, since this term was already employed in the calculation of the decay 
widths in potential model. In general the widths, which are the imaginary parts of the self- 
energies, can be a few hundred MeVs, which indicates that the chiral radiative corrections 
to the resonance masses cannot be small, and so should be included in computation of the 
masses. 

We shall now consider the corrections due to the chiral term (jSJ to the energy of an 
eigenstate of Hq. Let us denote the eigenenergy and normalized wavefunction by and 
\l'ni, respectively. Here m = {n,l,j,mj,q} denotes the set of quantum numbers classifying 
the eigenstate of the light quark, with n, q, and rrij denoting the radial excitation, quark 
flavor, and the angular momentum quantum numbers, respectively. The correction to the 
energy at one loop comes through the diagram in Fig. ^and is given by 



AE. 



where 



^ V- V- . f |jmn(fe)P 

and m^r denotes the mass of the light meson exchanged and represents the 5'f/(3)fiavor 
factors coming from the axial vertices. 

Using the Dirac equations for the wavefunctions the current can be written as 

Jmn(fc) = (ko -E^ + E^)p^^^{k) + p^^^ik) , (9) 
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FIG. 1: One loop correction to the energy of the eigenstate 



where 



P^^U^) = y rf'x ^L(5')7%5(m^ + m„ + 2K)^n(x)e*'^ (10) 
Substituting © into (|7j), and performing the integration over fco we obtain 

^E^ = ^^CnJ{m,^.m^) (11) 

n TT 

where 



.(m,n,..)^^/^^[(i.„-i..),p«^^^^ 
2Re[p^)(%^)(^)1. ^^^.^|^|^^';^_; 



(12) 



with E^{k) = + m2. 

We shall now write the currents pmn in terms of the radial functions of the eigenfunctions, 
which can be written as 

^n.(r)=| ^^"^-^'^ \y,,Ur), (13) 
gnijq{r)cr ■ r 

where yejmi'f') is the spinor harmonics. Since the light quark wavefunctions are eigenstates 
of the angular momentum operator, it is convenient to expand the plane wave exp{ik ■ r) in 
the definition of the currents pmn in (fTUI) as 

oo +e 

/k.r = ^^Y, ^'Jiikr) Yl Ye:Jf)YUk) ■ (W) 

^=0 m=-l 
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Then, the currents pmn"* can be expanded as 

pii'nH^) = T.P^^nm,i,m)YUk) (15) 
with, up to a common phase. 



= {jmjLm^\j'm'^)p'^l\\kll^) , (16) 

where 

Pmn = fmir)gn{r) - f^{r-)g^{r) , 

P^^l = {frn{r)gAr)+Ur)g^{r)){m^+m^ + 2Vs) . (17) 

Here {i,j,mj} are the angular momentum quantum numbers of the state m, {i',j',m'j} 
those of the state n and £^ is the angular momentum of the intermediate meson appearing 
in the expansion (IHj), and ■ r\\i'j') denotes the reduced matrix element. Eq. lfT^ 

provides a selection rule for the possible intermediate light meson angular momentum, 
for a given internal state and vice versa. 

Now, doing the angular part of the k- integration, which can be easily carried out with 
the decomposition (fTH|) . and using the unitarity relation 



we can rewrite the loop corrections to the energy as 



n ndTT 



where 

-.2 r u2 



J(m,n,g = y ^^^p„-i?m)|p«(|A:|,U| 

+2Re[p« (|^M.)p(^l*(|fcM.)] + ^ . 

We now focus on the energy corrections for the lowest energy parity doubling states, 

D(0-),Z),(0-) and D(0+),Z^,(0+). 



(20) 



They have the quantum numbers m = {1,0, |,±i,g = {d,s)}, {l,l,|,±|,g = {d,s)}, 
and, to shorten the notations, will be denoted by O^^s, 1^^^, respectively. For these states 
Pm.nd^l^^Tr) in ((201) are given by 

/•oo 

pS„(I^I,4) = / r'dr{Ur)g^ir)-Ur)g^ir))j,^{kr), (21) 

^0 

POO 

pSm(l^l,4) = / r2rfr(/^(r)^?„(r) + /„(r)^?^(r))(m^ + m„ + 2K)j,,(A;r).(22) 

Jo 

Before giving the numerical result, we comment on the divergence of the loop corrections. 
The loop correction J(m, n, for given m, n is free from ultraviolet (UV) divergence, with 
the wavefunctions providing the UV cutoff. However, the total loop correction obtained by 
summing over the internal states is quadratically divergent. The quadratic divergence comes 
from the first two terms of ()12|) . which can be easily summed over the internal states using 
the Dirac equation for the wavefunctions and the definition of pmn in (llUp. This gives the 
sum of the first two terms as 

^^cjuad.div. ^ ^^^Z^ f d'k j d^^^>lf{m^ + m^ + 2V,)^>^. (23) 

^ ^Jn [J {2'KfE^{k)\ J 

The third term in (jl2j) is at most linearly divergent. This quadratic divergence of the energy 
correction is not unexpected since the chiral quark model is an effective theory valid only 
at low energies. To regularize the UV divergence we introduce a three-momentum cutoff of 
the form 

e-^VA?,v (24) 

to the integrand in (pUj) . We regard Auv as the physical cutoff of the chiral quark model, 
but we shall see that our main result on the mass Gap depends little on the cutoff. 

To obtain the eigenfunctions of Hq we should first fix the parameters of the model. In 
the following we shall follow the setup as well as use the parameter values given in Ref . , 
in which the vector potential is Coulombic and the scalar potential is a linearly confining 
potential. The model has nine free parameters that are to be fixed by a global fitting of 
the predicted resonance masses to those observed values. For details we refer the readers 
to the above reference. Of course, the parameters were fixed without taking the chiral 
corrections into account, but we can use those values to estimate the loop correction effects 
at leading order, which are comparable in magnitude to the 1/Mc corrections and so play 



only a subleading role in fitting the parameters. We also need to fix the axial coupling gA, 
which is a free parameter in the chiral quark model. In principle it can be determined by 
fitting the hadronic decay width of an excited heavy-light meson to experimental data, or 
by lattice simulation. These approaches estimate qa to be around unity with a considerable 
uncertainty 7]. 

It is convenient to organize the energy corrections in terms of the angular momentum of 
the intermediate light mesons. For a given we sum over all allowed internal states, which 
can be selected by the angular momentum and parity conservations at the axial vertex, up 
to the first 10 radial excitations. As can be seen in Table H] the corrections drop rapidly at 
higher radial excitations. Our result is summarized in Table UTl at varying cutoffs Ayy. At 
a smaller cutoff Auv = 700 MeV the corrections drop quickly as increases, whereas at a 
larger cutoff Auv = 1200 MeV they drop slowly, reflecting the UV divergence of the chiral 
corrections. At all the cutoffs considered, the largest contributions come from the 1^^ = 1 
modes, and at larger cutoffs contributions from l.^ as large as four are signiflcant. 

Not surprisingly, the total energy corrections depend strongly on the UV cutoff. At 
Auv = 700 MeV they are a few hundred MeVs but at Auv = 1200 they are in GeV order. 
This shows that in our model the physical cutoff should be about 700 MeV. Although the 
total corrections are sensitive on the cutoff, we expect the difference of the interparity mass 
gaps between the parity doubling states is less sensitive on the cutoff. Indeed, summing the 
contributions up to !■„ = 9 we flnd that 

AEld-AEOd = -146,-312,-447 MeV, 

AEl,-AEOs = -271,-442,-579 MeV, (25) 

and 

Gap = {AEld - AEOd) - (AEl, - AEOs) 

= 125,130,132 MeV (26) 

at Auv = 700, 1000, 1200 MeV, respectively, and qa = 1- This shows that the chiral correc- 
tions shrink the interparity gaps, both in strange and nonstrange systems, but do so more in 
the strange system. This may be an explanation for the unusually smaller mass of Ds(2317) 
than given in the potential model. More interestingly, the mass Gap is remarkably stable 
under variation of the cutoff. This suggests that the Gap in Eq. comes almost entirely 

8 



1 2 3 4 5 6 

k(fm"') 



FIG. 2: The differential Gap G{k). No cutoff applied (Auv = cxd) and qa = 1- Notice that the 
bulk of the contribution to the Gap comes from the low energy region. 

from the low energy region far down the cutoff. To see this we plot in Fig. 2 the differential 
Gap Q{k), defined as the Gap before the integration over the momentum variable k, that is, 

POO 

Gap = / Q{k)dk , 

^0 

which can be obtained from the proper combination of the integrands for the energy cor- 
rections given in Eq. (j^IH) . The discontinuity in the plot comes from our implementation 
of the principal value prescription for the last term in the integrand in Eq. ()2()j) which 
has a pole at kp = ■>/ (-Em ~ E^)'^ — when the external state is 1^ and internal state is 
Od (^m = 2282 MeV, = 1895MeV,m^ = 140 MeV). Our numerical code handles the 
principal value integration using the identity 

Jo At Atp Jo J2kp "'P 

which is valid for any smooth function f{k). With this the pole at kp is now removed from 
the integrand and there appears a discontinuity at 2kp. Notice that the differential Gap has 
a peak around k ^ 250 MeV and the bulk of the contribution to the Gap comes from the 
low energy region. This shows that our prediction of the Gap is not affected by the UV 
physics and thus safe from the truncation of the higher derivative terms in the chiral quark 
model. 

Without the chiral corrections the potential model predicts almost vanishing Gap, while 
experimentally it is about 95 MeV. If we take = 0.82, as given in Ref. 0, we get 
experimentally consistent 90 MeV for the Gap. We note that a similar result was observed 



TABLE I: Energy corrections from the first five radial excitations at = 0. Gap = (AEld — 
AEOd) - {AEls - AEOs). Values are at Auv = 700 MeV and qa = 1- Units are in MeV. 



n 


1 


2 


3 


4 


5 


AEOa 


-23 


-3 











AEU 


29 


-49 


-8 


-1 





AEOs 


-33 


-3 











AEls 


-108 


-49 


-9 


-1 





Gap 


127 





1 









by Eichten using the heavy-light chiral Lagrangian but this was obtained without taking 
into account the tree-level mass terms for the heavy-light mesons arising from the explicit 
chiral symmetry breaking. 

It is notable that the Gap is dominated by the contributions from the n = 1, = 
modes, as can be seen from the Tables II|TT1 Interestingly, these modes widen the interparity 
mass gap in nonstrange system whereas in strange system narrow the gap, but still the Gap 
from these modes are already consistent with experiment. 

In conclusion, we calculated the one-loop chiral corrections for the heavy-light mesons in 
potential model based on the truncated chiral quark model, and have shown that the chiral 
corrections can account for the unusually small mass of Ds(2317) and the narrow mass 
difference between Ds{2317) and D(2308). Our calculation strongly supports the two-quark 
picture of the new resonances composed of a heavy quark and a light valence quark. 
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this work has been done. 



[1] BABAR Collaboration, B. Aubert et. ai, Observation of a narrow meson decaying to dfn* 
at a mass of 2.32-gev/c^, Phys. Rev. Lett. 90 (2003) 242001. 

[2] CLEO Collaboration, D. Besson et. al., Observation of a narrow resonance of mass 

2.46-gev/c**2 decaying to d/s*+ piO and confirmation of the d/sj*(2317) state, Phys. Rev. 
D68 (2003) 032002. 

[3] K. Abe et. al., Measurements of the d/sj resonance properties, Phys. Rev. Lett. 92 (2004) 
012002. 

[4] M. A. Nowak, M. Rho, and I. Zahed, Chiral effective action with heavy quark symmetry. 



11 



Phys. Rev. D48 (1993) 4370-4374. 
[5] W. A. Bardeen and C. T. Hill, Chiral dynamics and heavy quark symmetry in a solvable toy 

field theoretic model, Phys. Rev. D49 (1994) 409-425. 
[6] W. A. Bardeen, E. J. Eichten, and C. T. Hill, Chiral multiplets of heavy-light mesons, Phys. 

Rev. D68 (2003) 054024. 
[7] M. Di Pierro and E. Eichten, Excited heavy-light systems and hadronic transitions, Phys. 

Rev. D64 (2001) 114004. 
[8] Belle Collaboration, K. Abe et. al. Study of h- d**0 pi- (d**0 d(*)+ pi-) decays, 

Phys. Rev. D69 (2004) 112002. 
[9] H.-Y. Cheng and W.-S. Hou, B decays as spectroscope for charmed four-quark states, Phys. 

Lett. B566 (2003) 193-200. 
[10] K. Terasaki, Bahar resonance as a new window of hadron physics, Phys. Rev. D68 (2003) 

011501. 

[11] T. Barnes, F. E. Close, and H. J. Lipkin, Implications of a d k molecule at 2.32-gev, Phys. 

Rev. D68 (2003) 054006. 
[12] E. van Beveren and G. Rupp, Observed d/s(2317) and tentative d(2030) as the charmed 

cousins of the light scalar nonet, Phys. Rev. Lett. 91 (2003) 012003. 
[13] R. N. Calm and J. D. Jackson, Spin-orbit and tensor forces in heavy-quark light-quark 

mesons: Implications of the new d/s state at 2.32-gev, Phys. Rev. D68 (2003) 037502. 
[14] J. L. Goity and W. Roberts, A Relativistic chiral quark model for pseudoscalar emission 

from heavy mesons, Phys. Rev. D60 (1999) 034001. 
[15] A. Manohar and H. Georgi, Chiral quarks and the nonrelativistic quark model, Nucl. Phys. 

B234 (1984) 189. 

[16] E. Eichten, Talk presented at International workshop on heavy quarkonium, Beijing, Oct. 
12-15, 2004. 



12 



